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ABSTRACT 



We report on a detailed radiative transfer modeling of the observed scatter- 
ing polarization in the Ha line, which allows us to infer quantitative informa- 
tion on the magnetization of the quiet solar chromosphere. Our analysis sug- 
gests the presence of a magnetic complexity zone with a mean field strength 
(B) > 30 G lying just below the sudden transition region to the coronal tempera- 
tures. The chromospheric plasma directly underneath is very weakly magnetized, 
with (B) ~ 1 G. The possible existence of this abrupt change in the degree of mag- 
netization of the upper chromosphere of the quiet Sun might have large signifi- 
cance for our understanding of chromospheric (and, therefore, coronal) heating. 

Subject headings: magnetic fields — polarization — radiative transfer — scatter- 
ing — Sun: chromosphere — Sun: corona 



1. Introduction 



A long-standing issue in solar astrophysics concerns the strength and structu re variations 



with height o f the magnetic field in the chromosphere of the quiet Sun (e.g., IJudgd 12006 



Harveyl 120091 ) . Our empirical knowledge on this issue has remained vague notwithstanding 



the qualitative information provided by high resolution monochromatic images of the solar 
atmosphere taken at various wavelengths across strong spectral lines like Ha. For example, 
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Ha line-center images show a mass of c ell-spanning f ibrils as a flattened carpet, with upright 
ones jutting out from network patches (jRuttenll2007l ). Unfortunately, such images at various 
wavelengths (A) across the intensity profiles (/(A)) of chromospheric lines do not provide 
quantitative information on the magnetic field vector. 

The t emperature mini mum region of solar atmospheric models is transparent to Ha 
radiation (jSchoolmanl Il972l ) . As a result, we see the chromosphere at the very line center 
of Ha but the photosphere in the line wings. It is thus not surprising that the response 
function of the emergent circular polarization induced b y the Zeeman effect to magnetic fiel d 
perturbations shows large photospheric contributions (jSocas-Navarro fc Uitenbroek! 12004 ) . 
On the contrary, in the quiet Sun the linear polarization observed in Ha and in many 
other spectral lines is fully dominated by the presence of radiatively induced population 
imbalances and qua ntum coherences among the magnetic sublevels of the line's levels (e.g., 
Trujillo Buenol 120091). which produce Q/I profiles whose ma ximum values are usually located 
at the line center (jStenflo fc Kellerl 119971 ; iGandorferl 120001 ) . Moreover, this scattering line 
polarization is modified by the Hanle effect, which operates mainly in the line core and gives 
rise to Q/I and U/I profiles different to those cor responding to the zero-field case (e.g., 
Stenflolll994l ; lLandi DegFInnocenti fc Landolfil 120041 ) . Therefore, the Hanle effect in strong 
lines like Ha is the physical mechanism that should be exploited for facilitating quantitative 
explorations of the magnetism of the "quiet" solar chromosphere. 

In the quiet Sun the fractional linear polarization signals of the Ha line are weak (i.e., 
of the order of 0.1% when observing close to the solar limb), so that their detection with 
the present telescopes requires sacrifice the spatio-temporal resolution to be able to achieve 
high polarimetric sensitivity. As a result, it is easier to detect Stokes Q (with the tangent 
to the nearest solar limb being its reference direction) than Stokes U signals (whose positive 
and negative values quantify the rotation of the direction of linear polarization). For these 
reasons, typically only the St o kes Q /I profile is available, as is indeed the case with the 
Ha observation by IGandorferl (120001 ) that will be interpreted here. One could thus think 
that without U/ 1 it is presently impossible to infer the presence of an unresolved magnetic 
field in the quiet solar chromosphere, because it would require confronting the observed Q/I 
amplitude in the Ha line with the one that the highly inhomogeneous and dynamic solar 
chromospheric plasma would produce if it were unmagnetized. This strategy could be applied 
successfully for determining the mean magnetization of the quiet solar photosphere by solving 
the radiative transfer problem for the Sri at 4607 A line i n a realistic three-dimension al (3D) 
hydrodynamical model of the quiet solar photosphere (jTrujillo Bueno et al.ll2004l ). but a 
similar approach for inferring instead the magnetization of the quiet chromospheric plasma 
is not yet possible mainly because it is still computationally prohibitive to produce a realistic 
3D model of the thermal, density and dynamic structure of the quiet chromosphere (e.g., 
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Leenaartal2010l ). Fortunately, the Q/I profile observed by lGandorferl (120001 ) in a quiet region 
at about 5" from the solar limb (i.e., at fi = cos 0~ 0.1, with 9 being the heliocentric angle) 
shows a peculiar line core asymmetry (LCA; see the wiggly line in any of the right panels of 
Fig. [3]) whose very existence cannot be explained by the mere fact that the Ha line results 
from the superposition of seven components, four of which contribute to Stokes Q (see Fig. [p. 
Moreover, the Q/I profile of the (photo-ionization dominated) Ha line is not very sensitive 
to the chromospheric thermal structure, which justifies the strategy adopted here of using 
various (hot and cool) one-dimensional semi-empirical models of the quiet chromosphere for 
investigating the physical ori gin of the observ ed LCA. It is also necessary to note that the 
observed intensity profile (see lGandorferll2000f ) . which is fully insensitive to the magnetic field 
of the quiet solar chromosphere, does not show any noteworthy line-core asymmetry and that 
the small atomic weight of hydrogen does not favor broadening by macroscopic fluid motions. 
As shown in this letter, the LCA present in the observed Q/I profile can be explained by 
the Hanle effect of a magnetic field in the solar atmosphere whose height variation suggests 
the presence of a magnetic complexity zone located in the upper chromosphere of the quiet 
Sun. 



2. The Radiative Transfer Problem 

In the absence of collisional depolarization the critical value of the magnetic strength (in 
gauss) for the onset of the Hanle effect in a given atomic level is Bh=1- 137x 10~ 7 / (tnf e g), with 
g the level's Lande factor and tnf e its radiative lifetime in seconds (see the legend of Fig.[T|). 
Typically, the Hanle effect is sensitive to magnetic strengths between 0.1B H and 10B H , but 
note that due to collisional quenching the real Bh values of the nlj levels of hydrogen can be 
significantly larger. In fact, collisions with protons and electrons are efficient in producing 
transitions between the (n, and the (n, Z±l, j ) levels. Such collisional transitions play a 
non-negligible depolarizing role because at some heights their rates are similar to the inverse 
radiative lifetimes of the hydrogen levels. For B > 10 B H (j u ) (with B H (j u ) being the B H 
value of the line's upper level) the Hanle effect is only sensitive to the orientation of the 
magnetic field vector, but not to its strength (the saturation regime of the Hanle effect). 

Since the quiet solar chromosphere is optically thick at the center of the Ha line, we need 
to take into account radiative transfer (RT) effects. To compute the emerge nt Q/I profile 



in models of the quiet solar atmosphere we have applied a computer program (jStepanll2008l ) 
that takes into account the overlapping of the hydrogen transitions when solving jointly the 
statistical equilibrium equations for the multipolar compon ents of the atomic density matrix 



and the Stokes- vector transfer equations (see Chapter 7 in lLandi Degl'Innocenti fc Landolfi 
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2004J) via a fast iterative method and accurate formal solver (iTrujillo Buenol 120031 ). As 
justified in §1, we have used various (hot and cool) one-dimensional semi-empirical models 
of the quiet solar atmosphere assuming that at each height, z, there is a random-azimuth 
magnetic field of strength B(z), inclined by an angle 9b(z) with respect to the local vertical. 
Our conclusions on the physical origin of the observed LCA are based on multilevel RT 
calculations in such models, and it suffices with s howing in the followi ng section the results 
for the well-known FAL-C semi-empirical model ( iFontenla et al.l 119931 ). The chosen atomic 
model includes all the fine structure levels corresponding to the first four n-levels of the 
hydrogen atom. The neglect of hyperfine structure (HF S) should not produce any sizabl e 
overestimation of the scattering polarization amplitudes (jBommier fc Sahal-Brechotlll982l ). 
We quantify the excitation state of each nlj level by means of the multipolar comp onents 
of its corresponding atomic density matrix (ILandi Degl'Innocenti fc Landolfil 120041 ) . Such 
Pn(nlj) elements (with < K < 2j and —K < Q < K) quantify the level's overall 
population the population imbalances among its sublevels (pq >0 ) and the quantum 
coherences between each pair of them (Pq^q). We neglect quantum coherences between 
different j levels. This is a good approximation for investigating the Hanle effect in the Ha 
line, assuming that the Hanle effect is produced by the only action of magnetic fields of 
strength 5<100G. The reason is that under such circumstances the selection rule A/ = ±1 
for radiative transitions inhibits radiative couplings between levels (n,l,j) and (n, Z±l,j'). 
Moreover, for 5<100 G the energy separation between levels (n, l,j) and (n, l±2,j') is larger 
than both their natural width and than the ensuing Zeeman splitting. 

As mentioned above, the consideration of isotropic collisions with protons and electrons 
is important for a rigorous modeling of the scattering polari zation in hydrogen lines. We 
have applied the semi-classical impact approximation theory (jSahal-Brechot et al.lll996l ) for 
calculating the dominant dipolar collisional rates between the fine-structure levels. The 
Stark broadenin g of the hydr ogen lines is accounted for by using a suitable approach for the 
hydrogen lines (jStehlel 119961 ) . The inelastic collisiona l rates between the different n-levels 
were also taken into account (IPrzybilla fc Butlerll2004al ). With these physical ingredients our 
non-LTE synthesis of the Ha intensity profile is in good agreement with that obtained by 



other researchers (e.g.. IPrzybilla fc Butlerll2004bl ) using an atomic model with more n-levels 



Finally, it is also important to clarify that our calculations a re based on the complete fre- 
quen cy redistribution (CRD) theory of spectral line polarization (ILandi Degl'Innocenti fc Landolfi 
20041 ) which might be considered unsuitable for the Ha line given that the electron and pro- 
ton densities are sensitive to partial redistribution effect s through multilevel in terlockings 
with the other transitions included in our atomic model (IHubeny fc Liteslll995l ). However, 
in our calculations we have used the electron and proton densities of the FAL-C model itself 
and the Lyman lines tend to be optically thick in the atmospheric region of the Ha line for- 
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mation. Under suc h circumstanc es the CRD approximation is a suitable one for subordinate 



lines like Ha (e.g., iHeinzell Il995h 



3. The Physical Interpretation 

The emergent fractional linear polarization Q(X)/I(X) profile of the Ha line is basically 
the result of the superposition of four Qij(X)/I(X) blended components, where /(A) is the 
total intensity profile and each Qij(X) profile can be obtained by solving the transfer equation 
for Stokes Q using the /-component of the total absorption coefficient (rji) and the compo- 
nent's emissivity 6q (which is symmetric with respect to the component's central wavelength): 
Qij(X)~ J oo (eg/ry/)e~ t dt, with t the total optical thickness along the line of sight (LOS). As 
the fine-structure components i—j of the Ha line are shifted from line center (see Fig.[T]), the 
optical depth scale is not symmetric with respect to the component central wavelength A^- 
(i.e., r)i(Xij — AA) 7^ T]i(Xij + AA)). It follows that, in general, Qij(Xij — AA) ^ Qij(Xij + AA) 
so that the individual linear polarization components of the Ha line are asymmetric with 
respect to Ay. In conclusion: (1) The Qij/I profiles are asymmetric even in the absence of 
magnetic field; (2) the formation depth of the Qij/I components is a non-trivial function of 
wavelength. 

Figured shows that within the framework of our radiative transfer modeling the LCA 
in the observed Q/I profile cannot be explained neither in the absence of magnetic fields 
nor in the presence of homogeneous magnetic fields. The reason is that the shapes of the 
four Qij/I components are smooth "skewed gaussians" whose superposition gives rise to a 
rather smooth Q/I profile without any sizable asymmetry. The above-mentioned second 
conclusion provides the clue for our explanation of the observed LCA. In agreement with 
our numerical experiments, the LCA can be created by a modification of the Qij/I profiles 
of some components (mainly the "1" and "2" transitions of Fig.[2]) induced by a spatially 
varying magnetic field. The superposition of the blended linear polarization profiles can then 
give rise to a sizable LCA whose shape depends sensitively on the spatial distribution of the 
chromospheric magnetic field. 

Given that magnetic field models with a strength that decreases with height are also 
unable to explain the observed LCA we start with the following height-variation of the 
magnetic strength: 

B(z) = B ° w + B b . (1) 

v ' i _|_ e -{z-z )/w 10 v y 

This function, whose first term is a sigmoid, tends to B Q + Bb for z — > +oo (upper chromo- 
sphere) and to Bb for z — > — oo (lower photosphere). Note that B(zq) = Bb + B /2 and that 
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the smaller the positive parameter w the steeper the slope of the sigmoid. 

We fix for simplicity the inclination of our random-azimuth magnetic field and set B = 
50 G because in order to model the observed Q/I amplitude we need to achieve for the Ha 
line practically full saturation of the Hanle effect in the upper chromosphere. From non- 
LTE model calculations in various semi-empirical models we have found that zq and w have 
to be adjusted such that the field intensity in the upper chromosphere increases abruptly 
with height. For instance, in the FAL-C model the magnetic strength must increase from 
3 G at 1900 km to 50 G at 2200 km, approximately. Note that in this semi-empirical model 
the transition region to the coronal temperatures starts at a height z ~ 2200 km and that 
magnetic fields at larger heights in this model do not affect the linear polarization of the 
Ha emergent radiation for a LOS with \i = 0.1 because this line is already optically thin 
there. The other parameters of the magnetic field model are 9b and Bf,. Even though a 
horizontal field with a strength that increases with height in the upper chromosphere can 
induce a LCA, it does not seem to produce an optimal fit. We found a good agreement with 
the observed profile for 9b = 65°, zq = 2100km, w = 40km and _B;,~3G (see Fig.[3]A.). A 
comparison of the resulting Q/I profile with the observed one can be found in Fig.[3B- 

Although the agreement between the calculated and observed Q/I profiles is very good 
around the line cente r, there is now disagr e ement in the line wings, which are formed deeper 



in the atmosphere. iTrujillo Bueno et al.l (120041 ) showed that the bulk of the quiet solar 



photosphere is teeming with a distribution of tangled magnetic fields having a mean field 
strength (B) « 60 G (when estimating (B) assuming the simplest case, adopted here, of a 
single value field strength). We have included this important physical ingredient by adding an 
extra sigmoid function in Eq. ([I]) with a positive sign in the exponent, in order to get the B(z) 
variation given by Model 1 of Fig.[3]C. This model is still characterized by a magnetic zone 
with (B) > 30 G in the upper solar chromosphere, but has a strongly magnetized photosphere 
and a weakly magnetized lower chromosphere. The ensuing Q/I profile is shown in Fig.[3p. 
As can be seen, the consideration of this additional ingredient to satisfy the constraint of a 
"hidden" magnetic field in the quiet solar photosphere leads to a slightly better overall fit 
to the fractional linear polarization observed in the Ha line. We point out that this model 
for the height variation of the magnetic strength in the quiet solar atmosphere is not the 
only one that leads to a very good fit to the LCA of the observed Q/I profile. The key 
point is however that all plausible models are qualitatively similar, in the sense that all 
of them are characterized by the presence of a significantly magnetized zone in the upper 
solar chromosphere (where (B) must be at least 30 G), overlying weakly magnetized regions 
with (B) ~ 1 G. Model 2 of Fig.[3C shows a qualitatively similar magnetic field variation 
that leads almost to the largest possible further improvement in the Ha line wings without 
deteriorating the LCA fit (see the ensuing Q/I profile in Fig.|3p). 
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4. Concluding Comments 

The amplitude and shape of the scattering polarization profile of the Ha line is very sen- 
titive to the strength and structure of the magnetic field in the upper solar chromosphere. 
For this reason, we think that full Stokes-vector observations of the Ha line within and 
outside coronal holes would probably lead to an important breakthrough in our empirical 
understanding of chromospheric magnetism. A suitable instrument with which we are car- 
rying out such observations (in close collaboration with Dr. M. Bianda and Dr. R. Ramelli, 
from IRSOL) is the Zurich Imaging Polarimeter (ZIMPOL) attached to the Gregory Coude 
Telescope of the Istituto Ricerche Solari Locarno (IRSOL) and to the telescope THEMIS of 
the Observatorio del Teide (Tenerife; Spain). 

Our qne-dimensiqnal r adiative transfer interpretation of the LCA in the Q/I profile ob- 



served by iGandorferl (120001 ) suggests the presence of an abrupt magnetization in the upper 
chromosphere of the quiet Sun. In all the semi-empirical models of the quiet Sun atmosphere 
we have considered, the base of this magnetic complexity zone needed to explain the observed 
LCA lies a few hundred kilometers below the height where the model's sudden transition 
region to the coronal temperatures is located. Although the solar chromosphere is highly in- 
homogeneous and dynamic, we believe that our suggestion should be taken seriously because 
Ha is a photoionization-dominated line whose observed Stokes /(A) profile (which is fully 
insensitive to the chromospheric magnetic field) does not show any noteworthy asymmetry 
within the line core interval where the observed Q/I profile shows the peculiar asymmetry 
that motivated this investigation. 

The probable existence of an abrupt change in the mean magnetization of the upper 
chromosphere of the quiet Sun may have large significance for the passage of e merged magne- 



tized plasma across the chromosphere and into the corona (e.g., the review by iMoreno-Inseris 



20071 ) and for the overall energy balance of the outer solar atmosphere. For instance, coro- 
nal reconnection and dissipation of magnetic energy may be significantly modulated by the 
physical conditions in the upper chromosphere and by the h eight at which neighboring mag- 



netic fibrils begin to press against each other (jParkeri 120071 ). The presence of a strong and 
abrupt magnetization in the upper chromosphere of the quiet Sun might also hold the clue 
for clarifying the physical origin of the elusive acceleration mechanism of the solar wind. 
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3d 5/2 
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Fig. 1.— Components of the Ha line. Only 4 (solid lines) out of 7 (solid and dotted 
lines) blended transitions make a meaningful contribution to the scattering polarization 
of the Ha line in the solar atmosphere. The atomic polarization of the 2p%/2 level turns 
out to be negligibl e in the Ha formation region; hence selective absorption of polarization 
components (e.g., iManso Sainz fc Trujillo Buenol 120031 ) does not contribute to the linear 
polarization of Ha in the quiet Sun. The Einstein A u \ coefficients are given in s _1 . The arrows 
indicate the collisional dipolar depolarizing transitions. The three upper levels with angular 
momentum j > 1/2 have the following critical Hanle fields: -Bf/(3p 3 / 2 )=16 G, B H (3d 3 / 2 )=9 G 
and B H (3d$/2)=6 G, where B H is the magnetic strength for which the level's Zeeman splitting 
equals its natural width. 
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Fig. 2. — Q/I profiles of the Ha line and of its four polarizing components calculated in 
the FAL-C model atmosphere for a LOS with \i — 0.1. In each panel the full Q/I profile 
includes the contribution of the continuum polarization. The positive reference direction for 
Stokes Q is the parallel to the nearest limb. Left panel: the B = 0G case. We point out 
that identical Q/I shapes are obtained in the presence of a random-azimuth magnetic field 
with a strength that produces saturation of the upper- level Hanle effect (e.g., a horizontal 
field with B>50 G). The only difference is that the amplitude of the emergent polarization is 
smaller by a factor (3 cos 2 9b — l) 2 /4 (e.g., it is 0.12% for a horizontal magnetic field). Right 
panel: the emergent Q/I profiles in the presence of a random-azimuth magnetic field with 
constant strength B = 10 G and inclination 6b = 65°. Note that in this case the "2" and 
"3" transitions play now a more important role, in addition to transition "1" . In particular, 
transition "2" dominates the blue part of the full Q/I profile while transition "3" decreases 
the profile on the red side. As a result, the whole Q/I profile is slightly skewed toward 
blue wavelengths. However, no LCA similar to the observed one is produced because all 
components are depolarized in a rather wavelength-independent manner. 
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Fig. 3. — Magnetic field models (left panels) and calculated vs. observed Q/I profiles (right 
panels). The wiggly line in the right panels shows the Q/I profile of the Ha line observed 
by Gandorfer (2000) in a quiet region at about 5" from the solar limb. Panel B shows the 
calculated Q/I profile corresponding to the magnetic field model of panel A, while panel D 
shows the calculated Q/I profiles corresponding to the magnetic field models of panel C. 
All these magnetic field models produce a LCA similar to that present in the observed Q/I 
profile. Note in the left panels that for a LOS with \i — 0.1 the scattering polarization of the 
Ha line is insensitive to the structure of the magnetic field in the lower atmospheric region 
indicated by the dotted line part of the models. 



